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ABSTRACT 

A  comprehensive  computer  graphics  program  was  designed 
for  students  of  the  Naval  Postgraduate  School  to  assess  the 
vulnerability  of  air  targets  to  externally  detonating 
missile  warheads.  This  computer  simulation  program  was 
written  in  the  C  language  and  runs  on  the  IBM  PC/AT  or 
compatible.  It  illustrates  the  endgame  phase  of  a  typical 
air  defense  scenario  and  provides  the  students  the 
capability  to  comprehend  the  dynamic  situation. 
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I.  INTRODUCTION 

This  program  is  designed  to  simulate  the  endgame  phase 
between  a  missile  and  an  air  target  for  the  purpose  of 
illustrating  the  dynamics  of  the  encounter.  This  program, 
called  ENDGAME,  enables  the  user  to  determine  the  miss 
distance,  the  fragment  spray  zone,  and  the  best  point  for 
detonation  for  a  given  ordnance  package,  target  and 
missile  velocity  and   flight  path. 

This  program,  written  in  the  C  language,  and  using  the 
HALO  graphics  subroutine  packages,  provides  the  basic 
framework  from  which  further  simulations  of  increased 
complexity  and  sophistication  can  be  easily  implemented. 

To  establish  the  proper  perspective,  a  brief  outline  is 
given  as  follows;  the  following  chapter  gives  the 
theoretical  background  for  the  endgame.  This  is  followed 
in  Chapters  III  and  IV  by  a  construction  of  the  program  in 
the  C  language  with  comments  on  graphical  features. 
Finally  conclusions  are  drawn  in  Chapter  V,  and 
recommendations  are  made  for  future  improvements.  The 
program  variables  are  described  in  Appendix  A.  The  program 
listing  is  given  in  Appendix  B,  and  a  user's  manual  is 
given  in  Appendix  C. 


II.  THE  ENDGAME 


A.   THE  MISS  DISTANCE 

For  a  given  encounter  between  a  missile  and  a  target, 
the  guidance  system  can  use  any  one  of  several  methods  to 
navigate  the  missile  along  the  flight  path  from  launch  to 
the  intercept  with  the  target.  The  miss  distance  is 
basically  the  minimum  distance  between  the  missile  and  its 
target  during  the  encounter.  It  is  used  for  determining  the 
probability  that  a  fragment  from  the  detonation  of  the 
warhead  hits  the  target.  The  miss  distance  can  be 
evaluated  based  upon  the  encounter  conditions  shown  in 
Figure  2-1  and  assuming  the  velocity  of  both  the  missile 
and  the  target  are  constant. 
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Figure  2-1.   Encounter  Conditions  in  Grobal  Coordinates 
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To  evaluate  the  miss  distance,  grobal  coordinates  are 
chosen  for  convenience  and  proper  understanding.  See 
Figure  2-1. 

At  time  r 

T  +  t  =  M*+  m  +  s  (1.1) 

where 

T  =  the  vector  from  the  origin  to  the  target  at  r=0 

t  =  the  vector  from  the  target  at  r=0  to  the  target 
at  t=t 

-» 

M  =  the  vector  from  the  origin  to  the  missile  at  r=0 

-* 

m  =  the  vector  from  the  missile  at  t=0  to  the  missile 

at  t=t 
s  =  the  vector  from  m  to  t  (the  seperation  vector) 

Since  the  target  velocity,  Vt,  is  assumed  to  be  constant 

t  =  Vt  ti  (1.2) 

And  since  the  missile  velocity,  Vm,  is  also  assumed  to  be 
constant 


thus, 


m  =  Vm  (cosfl  l  +  sin^  3)  r  (1.3) 

s  =  (T  -  M)  +  (t  -  m) 

-  [(Tx  -  ly  +  Vt  r  -  Vm  COS^  r]  1         (1.4) 
+  [(Ty  -  My)  -  r  Vm  Sin^]  j" 

where  the  subscripts  x  and  y  denotes  the  x  and  y  components 


For  min  |s| 

d|s| 

=  0  (1.5) 

dr 

Solving  equation  (1.5)  for  r  using  equation  (1.4)  gives 
sx  <vm  cos*  "  vt)  +  sv  vm  sine 

T    =  ± (1.6) 

(Vm  cos*  -  Vt)2  +  (V  sin*)2 


Define         sx  =  Tx  ~  Mx 

Sy  =  Ty  -  My 
and  Vmtx  =  Vm  cos*  -  Vt 

vmty  =  vm  sin* 


Eq  (1.6)  can  be  written  as 


sx  vmtx  +  sy  vmty 
t  =  ± —  (1.7) 

vmtx  +  vmty 


and  eq  (1.4)  can  be  given  in  the  form 

2  =  [sx  -  vmtx  '3  1  +  Csy  "  vraty  '3  3         (1.8) 

Equation  (1.8)  can  be  used  to  determine  whether  the 
given  encounter  is  an  Early  Bird  or  a  Late  Bird  situation. 
In  the  Early  Bird  case,  the  missile  passes  in  front  of  the 
target,  whereas  it  passes  behind  the  target  in  the  Late 
Bird  case. 
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B.   THE  FRAGMENT  SPRAY  ZONE 

When  a  warhead  is  detonated  near  an  aircraft,  the 
fragments  are  usually  ejected  uniformly  and  begin  to 
propagate  outward  in  a  divergent  spherical  pattern  at  a 
velocity  that  is  the  vector  sum  of  the  initial  fragment 
velocity  from  a  static  warhead  detonation  and  the  missile 
velocity. [Ref .1] 

If  the  fragments  are  assumed  to  have  uniform  velocity 
and  to  be  uniformly  spread  over  a  spherical  segment,  the 
fragment  spray  density  at  the  distance  R  is  given  by 

N 

P = (1.9) 

27rR^(cos  ^  -  cos  <f>2) 

where 

N  =  the  total  number  of  fragments  in  the  warhead 

The  equation  defining  the  angles  <j>-^   and  <t>2   is 


<t>^   =  tan-1 


V_  sinfl  +  V  sin(  9   +  a*) 

10) 


Vm  COS0  +  VQ  cos ( 


{      0      f     Cti)  ,-J 

= -9  (1. 

+  a±)     ±   Vt  J 


i  =  1,2 
where 

Vm  =  the  speed  of  the  missile 
d     =  the  elevation  angle  of  the  missile 
VQ  =  the  average  fragment  speed  with  respect  to  a 
stationary  warhead 
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a^  =  the  static  leading  fragment  spray  angle 
a2  =  the  static  trailing  fragment  spray  angle 
V^  =  the  target  horizontal  speed 

C.   THE  BEST  POINT  FOR  DETONATION 

The  fragment  spray  density  from  the  detonation  point  is 
determined  by  the  given  detonation  distance.  At  this 
point,  the  following  problem  is  encountered.  Where  is  the 
best  point  for  detonation  ?  The  specific  answer  to  that 
question  requires  an  understanding  that  there  are  four 
different  fragment  spray  zones  created  by  the  encounter 
conditions  and  the  warhead  parameters  as  shown  in  Figure 
2-2. 


Figure  2-2.    The  Fragment  spray  zone 
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In  region  1  and  3,  part  of  the  fragment  spray  zone  hits 
part  of  the  target.  In  region  2,  the  entire  fragment  spray- 
zone  hits  part  of  the  target,  and  in  region  4,  part  of  the 
fragment  spray  zone  hits  all  of  the  target. 
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III.  PROGRAMMING  "ENDGAME" 

This  chapter  contains  the  programming  description  of 
the  implementation  of  the  endgame  equations  using  the 
Lattice  C  compiler.  Even  though  several  compilers  are 
available  these  days,  the  Lattice  C  compiler  is  used  since 
it  is  a  portable  compiler  for  the  high  level  programming 
language  called  C. 

A.  PREPARING  SOURCE  FILE 

The  source  file  was  written  by  using  the  text  line 
editor  EDLIN  in  IBM-DOS  version  3.2  to  create,  change,  and 
display  the  source  file.  All  the  external  variables  are 
declared  outside  of  the  module,  including  other  standard 
source  files.  The  driver  program,  main(),  is  the  first 
routine  to  be  executed.  This  is  followed  next  by  several 
subroutine  programs,  sub-subroutines,  and  so  on.  Since  the 
program  is  so  massive  that  it  is  a  little  hard  to  read  the 
whole  program  thouroughly  at  a  glance,  each  of  the 
subroutine  functions  are  entitled.  See  the  program  listing 
in  Appendix  B  for  further  particulars. 

B.  COMPILING 

The  Lattice  C  compiler  version  3.00  [Ref.3]  is  used 
for  compiling  the  source  file  into  the  object  file.  The  LC 
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command  is  used  with  several  options:  -w  option  to  shut  off 
warning  messages  generated  for  return  statements,  and  -k2 
option  to  specify  the  generate  code  for  80286.  Another 
option  chosen  is  -s  for  naming  program  segments  to  provide 
more  availability  on  RAM  working  memory.  Otherwise  the 
storage  requirements  exceed  the  available  memory. 

This  was  the  most  difficult  part  of  the  author's  thesis 
work:  to  find  the  proper  options  that  would  work  with  the  P 
model  as  the  program  got  bigger  and  bigger.  Refer  to  the 
Lattice  C  compiler  manual  for  more  detailed  information. 

C.   LINKING 

When  linking  the  object  file  into  the  program  file,  the 
LINK  command  in  MS-DOS  version  3.0  was  used.  Since  most 
parts  of  the  program  call  the  graphics  subroutine  functions 
in  HALO  to  draw  the  graphics  on  the  screen,  a  special 
format  for  linking  the  program  with  the  HALO  library  is 
required.  With  program  model  (version  3.0)  supported  by 
the  full  available  working  memory,  the  following  command  is 
used: 

LINK  CP+ENDGAME+HAL0DVP3 , ENDGAME , , LCMP+LCP+HALOP3 
Further  details  on  each  step  is  discussed  in  Section  4 
of  the  HALO  manual  and  can  be  refered  to  by  looking  into 
the  LINK  command  in  the  MS-DOS  manual. 
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IV.   GRAPHICS  IN  THE  ENDGAME  PROGRAM 

As  mentioned  in  the  Introduction,  most  of  the  author's 
efforts  were  devoted  to  creating  screens  that  illustrate 
the  endgame  geometry  using  colorful  graphics  to  obtain 
desirable  results  from  the  given  warhead  and  encounter 
data  input.  In  ENDGAME,  the  useful  information  is 
presented  on  each  screen  with  high  resolution,  is  processed 
and  displayed  graphically  on  the  screen,  simulating  each 
endgame  phase  and  showing  what  is  happening. 

All  the  screens  were  created  using  HALO  in  the  graphics 
mode,  except  the  option  screen  which  shows  the  sequential 
questions  and  input  data.  The  HALO  library  used  in  this 
work,  version  2.26,  is  a  collection  of  high  performance 
subroutines  which  allow  the  application  programmer  to 
develope  sophisticated  computer  graphics  programs. 

ENDGAME  uses  world  coordinates  throughout  the  program. 
The  world  coordinates  are  mapped  to  device  coordinates  on  a 
specified  range  from  0  to  640  on  x-axis  and  from  0  to  350 
on  y-axis  so  as  to  locate  each  point  within  the  resolution 
of  the  EGA (Enhanced  Graphics  Adaptor)  display  device.  The 
upper  left  corner  of  the  screen  is  at  coordinates  (0,0)  and 
the  lower  right  corner  is  at  coordinates  (639,349)  for 
convenience.  To  draw  objects,  the  coordinates  and 
position  of  the  objects  have  to  be  scaled  by  the  aspect 
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ratio.  Aspect  is  expressed  as  a  ratio  of  width  to  height. 
The  value  1.439  was  emperically  chosen  to  make  the  range 
increment  on  the  horizontal  axis  be  the  same  as  that  on 
the  vertical  axis. 

Another  factor  considered  was  trans_f actor  [Ref.2].  As 
mentioned  in  the  user's  manual  in  Appendix  C,  some 
particular  situations  had  to  be  excluded  in  the  endgame. 
The  user  should  exclude  the  encounter  conditions  where  the 
difference  in  relative  position  between  the  target  and  the 
missile  is  larger  than  1,000  ft.  Hence,  the  difference 
always  lies  within  1,000  ft  and  the  graphic  resolution  is 
taken  account  into  by  trans_f  actor.  The  value  of 
trans_f actor  such  as  0.2,  0.3,  etc.,  is  emperically  chosen 
based  upon  the  value  that  gives  good  graphic  resolution 
on  the  screen. 

A.   CREATED  SCREENS 

All  the  screens  created  by  HALO  are  in  the  graphics 
mode.  There  are  three  different  screens  which  allow 
students  to  comprehend  the  dynamic  situation  between  the 
target  and  the  missile. 

1.    The  Opening  Screen 

When  ENDGAME  is  executed,  the  opening  screen  is 
displayed,  welcoming  the  student  to  the  endgame  assessment 
as  shown  in  Figure  4-1.  Hitting  any  key  allows  the  student 
to  proceed  to  the  first  step. 
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2.  The  Miss  Distance  Screen 

After  inputting  the  data  requested  by  the  option 
screen,  the  student  will  see  the  miss  distance  screen 
(Figure  4-2) .  This  screen  tells  the  student  what  the  miss 
distance  is  numerically  and  shows  how  it  looks 
vectorically  with  the  target  and  the  missile  vectors  shown 
on  the  screen.  It  also  prints  and  shows  whether  the 
encounter  situation  is  an  Early  Bird  or  a  Late  Bird.  Refer 
to  the  Appendix  A  for  the  other  variables  shown  on  the 
screen. 

3 .  The  Fragment  Spray  Zone  Screen 

This  screen  (Figure  4-3)  shows  the  fragment  spray 
zone  with  respect  to  the  target.  Unfortunately,  the  screen 
doesn't  show  the  different  fragment  vectors  due  to  leading 
and  trailing  fragment  spray  angle  with  different  colors. 
However,  the  reader  should  be  able  to  visuallize  the 
fragment  pattern.  Also  refer  to  the  Appendix  A  for  other 
variables. 

4.  The  Best  Point  for  Detonation  Screen 

The  student  will  see  where  the  best  point  for 
detonation  is  on  this  screen  (Figure  4-4) .  It  shows  the 
zones  which  are  appropriate  for  the  detonation  of  the 
warhead  for  the  given  encounter  situation  considering  the 
combination  of  the  dynamic  fragment  speed  with  respect  to 
the  target  and  the  missile  velocity  vector.  From  this 
screen,  the  student  can  go  back  to  see  the  previous  screens 
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simply  by  following  the  instructions  on  the  bottom  of  the 
screen.  Finally,  the  end  screen  will  appear  if  the  student 
follows  the  appropriate  instruction  on  the  bottom  of  the 
screen. 

B.   OPTION  SCREENS 

Several  questions  appear  on  the  screen  prompting  the 
student  to  enter  the  necessary  input  data.  This  data  is 
used  to  calculate  the  miss  distance  and  the  best  point  for 
detonation.  These  screens  are  not  in  the  graphics  mode,  but 
are  in  standard  text  mode.  Warning  or  error  messages  might 
appear  if  the  student  inputs  unreasonable  data.  The 
limitations  for  this  endgame  program  are  described  in  the 
user's  manual  which  is  given  in  Appendix  C. 
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Figure  4-1.   The  Opening  Screen 
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Figure  4-2.   The  Miss  Distance  Screen 
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Figure  4-3.   The  Fragment  Spray  Zone  Screen 
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Figure  4-4.   The  Best  Point  for  Detonation  Screen 
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V.   SUMMARY  AND  RECOMMENDATIONS 

This  computer  graphics  simulation  program  for  the 
endgame  between  a  missile  and  an  air  target  was  written  in 
the  C  language  using  the  Lattice  C  compiler  running  on  MS- 
DOS.  The  program  is  divided  into  three  major  categories: 
miss  distance,  fragment  spray  zone,  and  best  point  for 
detonation.  It  was  intended  to  visually  demonstrate  the 
dynamic  evaluation  for  the  terminal  events  in  the  encounter 
between  an  aircraft  and  a  high  explosive  proximity  fuzed 
warhead  on  the  missile. 

It  is  hoped,  finally,  that  students  who  are  more 
interested  in  this  field  will  continue •  to  develop  the 
present  ENDGAME  program  with  an  anxiety  and  an  eagerness. 
Especially,  a  more  realistic  endeavor  to  expand  the  variety 
of  options  within  the  program  to  account  for  any  encounter 
situation  is  strongly  recommended.  For  example,  adding  a 
target  acceleration  to  the  target  motion  would  allow  more 
realism.  From  a  programming  point  of  view,  a  solution  to 
the  memory  expansion  problem  is  also  desirable.  Hopefully, 
the  probability  of  kill  will  be  developed  in  future  study. 
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APPENDIX  A 


PROGRAM  VARIABLES 


The  variables  used  in  endgame  computer  program  are 
listed  as  follows; 


alphal  [degree] 
alpha2  [degree] 


static  leading  fragment  spray  angle 
with  respect  to  missile  axis 
static  trailing  fragment  spray  angle 
with  respect  to  missile  axis 


gammal  [degree] 
gamma2  [degree] 


angle  from  horizontal  to  dynamic 
leading  fragment  spray 
angle  from  horizontal  to  dynamic 
trailing  fragment  spray 


ht 


[degree] 


target  heading 


length  [ft] 


target  presented  length 


minsx   [ft] 
minsy   [ft] 
miss_distance  [ft] 
mx     [ft] 
my     [ft] 


x  component  of  miss  distance  vector 
y  component  of  miss  distance  vector 
miss  distance  or  closest  point 
x  component  of  missile  position 
y  component  of  missile  position 


num 


total  number  of  fragments  in  warhead 


phione  [degree] 
phitwo  [degree] 


dynamic  leading  fragment  spray  angle 
with  respect  to  missile  axis 
dynamic  trailing  fragment  spray  angle 
with  respect  to  missile  axis 


r_det   [ft] 


rho 


[frag/ ft-] 


detonation  distance 
fragment  spray  density 


sx 

sy 


[ft] 
[ft] 


x  component  of  difference  between 

missile  and  target  position 

y  component  of  difference  between 
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missile  and  target  position 


tau 
theata 
tx 
ty 

[sec] 
[degree] 
[ft] 
[ft] 

vftl 

[ft/ sec] 

vft2 

[ft/ sec] 

vftxl 
vftx2 
vftyl 
vfty2 

[ft/ sec] 
[ft/ sec] 
[ft/ sec] 
[ft/ sec] 

vm 
vmtx 

[ft/sec] 
[ft/sec] 

vmty 

[ft/ sec] 

vo 

[ft/sec] 

vt 

[ft/sec] 

minimum  time  for  miss  distance 
elevation  angle  of  missile 
x  component  of  target  position 
y  component  of  target  position 


dynamic  leading  fragment  speed  with 
respect  to  target 

dynamic  trailing  fragment  speed  with 
respect  to  target 


x  component  of  vftl 
x  component  of  vft2 
y  component  of  vftl 
y  component  of  vft2 


speed  of  missile 

x  component  of  missile  speed  with 

respect  to  target 

y  component  of  missile  speed  with 

respect  to  target 


average  fragment  speed  with  respect 
to  stationary  warhead 
target  horizontal  speed 
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APPENDIX  B 


PROGRAM  LISTINGS 


The  endgame  computer  simulating  program  written  in  C  is 
as  follows. 
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/******    define  the  all  external  variables   *******/ 


/******  include  another  source  file  ******/ 

#include  "stdio.h" 
# include  "stdlib.h" 
# include  "math.h" 

/**********  image  file  memorized  *********/ 

static  char  filel []={ "miss. pic" } ; 
static  char  file2 [ ]={"prob.pic"} ; 
static  char  file3 []={ "best. pic"} ; 


/**********  position  variables  ***********/ 

float  xl=0. 0,yl=0. 0,x2=640. 0,y2=350.0;  /*world  coordinate*/ 
float  gcx,gcy  ;  /*  graphic  cursor  position  */ 

float  gcxl , gcx2 , gey 1 , gcy2 ;   /*  for  varing  graphic  cursor  */ 
float  tcx,tcy;  /*  text  cursor  position  */ 

float  dxfdy;  /*  cursor  increment  */ 

float  bxl,byl,bx2,by2  ;  /*  box  creation  */ 

/***  variables  for  adjusting  to  screen  ***/ 

float  trans_factor  ; 

float  dist_f actor , dist_l , dist_2 ; 

/************  color  definition  ***********/ 

int  black  =0; 
int  blue  =  1  ; 

int  green  =  2  ; 
int  cyan  =  3  ; 

int  white  =  7  ; 

int  liblue  =  9  ; 

int  ligreen  =  10 

int  licyan   =  11 

int  lired    =  12 

int  lipurple  =  13 

int  liyellow  =  14 


/***********  system  variables  ************/ 

float  aspect=0.7249;         /*  ratio  of  width  to  height  */ 
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int  mode  =4;  /*  enhanced  graphics  adapter  mode  */ 

int  page=2 ;        /*  available  graphics  page  for  system  */ 

/************  drawing  variables  **********/ 

float  angl,ang2;  /*  starting  and  ending  angle  of  arc  */ 

int  width;  /*  line  width  */ 

int  style;  /*  line  style  */ 

float  radius;  /*  radius  of  arc  in  world  coordinates  */ 

char  s[15] ; 

int  c; 

int  h,w,p,m; 

float  itx,ity; 
float  imx,imy; 
int  count; 

/*********  variables  for  program  *********/ 

float  tx,ty  ; 

float  mx,my  ; 

float  sx,sy  ; 

float  vt,vm  ; 

float  vo  ; 

float  ht  ; 

float  theata; 

float  tau; 

float  vmtx,vmty  ; 

float  minsx,minsy  ; 

float  m_dist; 

float  vftxl,vftyl; 

float  vftx2,vfty2; 

float  gammal , gamma2 ; 

float  alphal,alpha2 ; 

float  vftl,vft2; 

float  phione,phitwo  ; 

float  r  det,num, length, rho; 
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/**  main  function  **/ 

/it********************************************************/ 


/*   This  states  that  it  is  the  first  or  main  routine  to  be 
executed.   It  shows  the  structure  of  the  endgame  program  */ 


main() 

{ 

setdev("haloibme.devH) ; 
setscreen(&page) ; 
initgraphics(&mode) ; 
setasp(&aspect) ; 
setworld ( &xl , &y 1 , &x2 , &y2 ) ; 


opening () ; 
optionl ( ) ; 
drawl ( ) ; 
option2 ( ) ; 
draw2 ( ) ; 
draw3 ( ) ; 


closegraphics ( ) ; 
} 


/*  define  graphics  system  */ 

/*  set  screen  type  */ 

/*  start  with  HALO  */ 

/*  set  aspect  ratio  */ 

/*  world  coordinates  */ 

/*  opening  screen  */ 

/*  miss  distance  screen  */ 


/*  fragment  spray  zone  screen  */ 
/*  best  point  for  detonation  screen  */ 

/*  end  up  with  HALO  */ 
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/it********************************.************************/ 

/**  SUBROUTINE:   OPENING ()  **/ 

/it********************************************************/ 

opening ( ) 

{ 

varl3 () ; 

varl4 ( ) ; 

setcolor(&lired) ; 
gcx=3  0.0, gcy=3  20.0; 
movabs ( &gcx , &gcy ) ; 
gcx=32  0.0; 
lnabs(&gcx, &gcy) ; 

settextclr (&lired, &black) ; 

tcx=3  4  0.0, tcy=3  20.0; 

movtcurabs (&tcx, &tcy) ; 

text ("Subject  :  AIR  DEFENCE  LETHALITY"); 

setcolor(Stlipurple)  ; 
gcx=4  0.0, gcy=3 10.0; 
movabs ( &gcx, &gcy) ; 
gcx=330. 0; 
lnabs(&gcx, &gcy) ; 

settextclr (&lipurple, &black)  ; 
tcx=3  50.0,tcy=310. 0; 
movtcurabs ( &tcx , &tcy ) ; 
text ("Professor  :  R.  E.  Ball") ; 

h=4,w=4 ; 

settext (&h, &w, &p, &m) ; 
settextclr ( &green, Sblack) ; 
tcx=2 10 . 0 , tcy=2  2  0.0; 
movtcurabs ( &tcx, &tcy) ; 
text ("WELCOME") ; 

h=3,w=3 ; 

settext ( &h , &w , &p , &m) ; 
settextclr ( Sligreen, Sblack)  ; 
tcx=299 . 0 , tcy=185 . 0 ; 
movtcurabs ( &tcx, &tcy)  ; 
text ("to") ; 

h=6,w=6; 

settext ( &h , &w , &p , &m) ; 
settextclr (&licyan, Sblack) ; 
tcx=155 . 0 , tcy=107 . 0 ; 
movtcurabs ( &tcx , &tcy ) ; 
text ("ENDGAME") ; 
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h=2,w=2 ; 

settext(&h, &w, &p, &m) ; 

settextclr(&liyellow, Sblack) ; 

tcx=123 . 0,tcy=55.0; 

movtcurabs(&tcx,&tcy) ; 

text ("NAVAL  POSTGRADUATE  SCHOOL"); 

bottom () ; 

varl2 () ; 

} 


/**  SUBROUTINE:  OPTIONl()  **/ 

/•A*******************************************************/ 

optionl() 

{ 

labl: 

printf("\n  Initial  x-position  of  target  [ft]  ?  :"); 

scanf ("%f",&tx) ; 

if((tx<0.0) | | (tx>100000.0) ) 

{ 

print f ("\n  Unreal!  Assume  that  it  be  above  ground  and 

real.")  ; 

printf("\n  Type  again."); 

goto  labl; 

} 
else 

{ 

lab2: 

printf("\n  Initial  y-position  of  target  [ft]  ?  :"); 

scanf ("%f",&ty) ; 

if((ty<0.0) | | (ty>100000.0) ) 

{ 

print f ("\n  Unreal!  Assume  that  it  be  above  ground  and 

real . " )  ; 

printf("\n  Type  again."); 

goto  lab2 ; 

} 
else 

{ 

lab3 : 

printf("\n  Target  speed  [ft/sec]  ?  :"); 

scanf ("%f",&vt) ; 

if((vt<0.0) | | (vt>2000)) 

{ 

printf("\n  Unreal!  Type  again."); 

goto  lab 3 ; 

} 
else 
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{ 

lab4: 

printf ("\n  Target  heading  [degree]  ?  (  type  0  or  180  ) 
■  it  \  • 

scanf ("%f", &ht) ; 

if ((ht!=0) &&(ht 1=180) ) 

{ 

printf ("\n  Assume  level  flight!  Type  again."); 
goto  lab4; 

} 
else 

{ 

lab5: 

printf ("\n  Initial  x-position  of  missile  [ft]  ?  :"); 

scanf ("%f",&mx) ; 

if ((mx<0.0) | | (mx>100000.0) ) 

{ 

printf ("\n  Unreal!  Assume  that  it  be  above  ground  and 

real  . ■ ) ; 

printf ("\n  Type  again."); 

goto  lab5 ; 

} 
else 

{ 

lab6: 

printf ("\n  Initial  y-position  of  missile  [ft]  ?  :"); 

scanf ("%f",&my) ; 

if((my<0.0) | | (my>100000.0) ) 

{ 

printf ("\n  Unreal!  Assume  that  it  be  above  ground  and 

real . " )  ; 

printf ("\n  Type  again."); 

goto  lab 6 ; 

} 

else  if (( (abs(tx-mx)<100.0)&&(abs(ty-my)<100.0) )  || 
((abs(tx-mx)>1000.0)&&(abs(ty-my)>1000.0) ) ) 

{ 

printf ("\n  Unproper  input  data  for  end  game!  Start 

again.") ; 

goto  labl; 

} 

else 

{ 

lab7: 

printf ("\n  Missile  speed  [ft/sec]  ?  :"); 

scanf ("%f",&vm) ; 

if ( (vm<vt) | | (vm>6000.0) ) 

{ 

printf ("\n  Unreal!  Missile  speed  must  be  greater  than 

target  speed."); 

printf ("\n  Make  it  real!  Type  again."); 
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goto  lab7; 

} 
else 

{ 
lab8: 

printf ("\n  Elevation  angle  of  missile  (theata)  [degree]  ? 

:"); 

scanf ("%f",&theata) ; 

if  ( 

( ( (ht==0) && ( ( ( (tx<mx) && (ty>my) ) && ( (180<=theata) && (3  60>=thea 

ta))) 

I  I 

(((tx<mx)&&(ty<my) ) &&( (0<=theata) && (theata<=180) ) ) ) ) 

I  I 

( (ht==0) && ( ( ( (tx>mx) && (ty>my) ) && ( (theata>=90) && (theata<=360 

))) 

| | ( ( (tx>mx) && (ty<my) ) && ( (theata>=0) && (theata<=270) ) ) ) ) ) 

I  I 

( ( (ht==180) && ( ( ( (tx<mx) && (ty>my) && ( (theata>=180) && (theata<= 

90))) 

I  I 

( ( (tx<mx) && (ty<my) ) && ( (theata>=270) && (theata<=180) ) ) ) ) 

II 

( (ht==180) && ( ( ( (tx>mx) && (ty>my) ) && ( (theata>=270) && (theata<= 

360))) 

I  I 

(((tx>mx)&&(ty<my) ) && ( (theata>=0) &&(theata<=180) ) ) ) ) ) 

)) 

{ 

printf ("\n  Target  might  not  be  hit.  Put  appropriate  angle 

i  I  n  \  . 

goto  lab8 ; 

} 
else 

{ 

ht  =  ht  *  PI  /  180.0; 
theata  =  theata  *  PI  /  180.0; 
sx  =  tx  -  mx  ; 
sy  =  ty  -  my  ; 
if (ht==0) 

vmtx  =  vm  *  cos (theata)  -  vt  ; 
else 

vmtx  =  vm  *  cos (theata)  +  vt  ; 
vmty  =  vm  *  sin (theata)  ; 
tau  =  (  sx  *  vmtx  +  sy  *  vmty  ) / ( 
pow (vmtx ,2.0) +pow (vmty ,2.0)) ; 
minsx  =  sx  -  tau  *  vmtx  ; 
minsy  =  sy  -  tau  *  vmty  ; 

m_dist  =  sqrt (pow (minsx ,2.0)  +  pow(minsy,2 . 0) ) ; 
ht  =  ht  *  180.0  /  PI  ; 
theata  -  theata  *  180.0  /  PI  ; 
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change : 

printf ("\n  Would  you  like  to  change  these  values  ?  (  Y/N  ) 

:")  ; 

c=getch() ; 

if(  c=='y') 

{ 

printf ("\n  Yes.  Hit  any  key  to  change."); 

varl2 ( ) ; 

goto  labl; 

} 
else  if(  c=='n') 

{ 
printf ("\n  Hit  any  key  to  continue."); 

} 
else 

{ 

printf ("\n  Type  again  correctly."); 
goto  change ; 

} 
var!2()  ; 

} 


/it********************************************************/ 

/**  SUBROUTINE:  DRAW1()  **/ 

/*********************************************************/ 

drawl ( ) 

{ 

top ( ) ; 

text ("MISS  DISTANCE  [ft]  =") ; 

inputl() ; 

inputll() ; 

coord_l ( ) ; 

loc_l ( ) ; 

vectorl ( ) ; 

indexl() ; 

output 1 ( ) ; 

outputll() ; 

output 12 ( ) ; 

} 
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/It********************************************************/ 

/**  SUBROUTINE:  OPTION2()  **/ 

/it********************************************************/ 

option2 () 

{ 

labll: 

printf("\n  Fragment  velocity  (Vo)  [ft/sec]  ?  :"); 

scanf ("%f",&vo) ; 

if  ((vo<0.0)  |  |  (VO200000.0)) 

{ 

printf("\n  Unreal!   Assume  that  it  be  real."); 

printf("\n  Type  again."); 

goto  labll; 

} 
else 

{ 

labl2: 

printf("\n  Spray  angle  (alphal)  [degree]  ?  s"); 

scanf ("%f",&alphal) ; 

if  ((alphaKO.O)  |  |  (alphal>180.  0) ) 

{ 

printf("\n  Unreal!   Type  again."); 

goto  labl2; 

} 
else 

{ 

labl3: 

printf("\n  Spray  angle  (alpha2)  [degree]  ?  :"); 

scanf ("%f"/&alpha2) ; 

if  ((alpha2<0.0) | | (alpha2>180. 0) ) 

{ 

printf("\n  Unreal!   Type  again."); 

goto  lab 13 ; 

} 
else 

{ 

labl4: 

printf("\n  Number  of  fragments  (N)  ?  :"); 

scanf  ("%f",&num)  ; 

if (num<0.0) 

{ 

printf("\n  Unreal!   Type  again."); 
goto  labl4 ; 

} 

else 

{ 

labl5: 

printf("\n  Detonation  distance  (R)  [ft]  ?  :"); 

scanf ("%f",&r_det) ; 

if  (r_det<0.0) 
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{ 

printf("\n  Unreal!   Type  again."); 
goto  labl5; 

} 
else 

{ 
labl6: 

printf("\n  Target  presented  length  (L)  [ft]  ?  :"); 

scanf ("%f",&length) ; 

if ((length<0.0) | | (llength>200.0) ) 

{ 

printf("\n  Unreal!   Type  again."); 

goto  labl6; 

} 
else 

{ 

alphal  =  alphal  *  PI  /  180.0  ; 
alpha2  =  alpha2  *  PI  /  180.0  ; 

varl5(); 

if (ht!=0.0) 

vt=-vt ; 
else 

vt=vt ; 
vftxl  =  viti  *  cos (theata)  +  vo  *  cos(theata  +  alphal)  - 

vftx2  =  vm  *  cos(theata)  +  vo  *  cos(theata  +  alpha2)  - 

vftyl  =  vm  *  sin(theata)  +  vo  *  sin(theata  +  alphal)  ; 

vfty2  =  vm  *  sin (theata)  +  vo  *  sin (theata  +  alpha2)  ; 

gammal  =  atan(vftyl/vftxl)  ; 

gamma2  =  atan(vfty2/vftx2)  ; 

vftl  =  sqrt(pow( vftxl, 2.0)  +  pow (vftyl, 2 . 0) ) ; 

vft2  =  sqrt (pow(vftx2,2 . 0)  +  pow(vfty2 , 2 . 0) ) ; 

varl6() ; 

phione  =  gammal  -  theata  ; 
phitwo  =  gamma 2  -  theata  ; 

rho=(num/(2*PI*(pow(r_det,2.0) ) * (cos (phione) - 
cos (phitwo) ) ) ) ; 

} 
} 
} 

} 
} 
} 
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vt 
vt 


Iabl7; 

print f ("\n  Would  you  like  to  change  these  values  ? 

(Y/N):»); 

c=getch() ; 

if   (C==Iyl) 

{ 

printf("\n  Yes.  Hit  any  key  to  continue."); 
varl2 () ; 
goto  labll; 

} 

else  if  (c=='n' ) ; 

{ 
printf("\n  Hit  any  key  to  continue."); 

} 
else 

{ 

printf("\n  Type  again  correctly."); 

goto  lab 17 ; 

} 

varl2 ( ) ; 

} 


/**  SUBROUTINE:  DRAW2()  **/ 

/*********************************************************/ 

draw2 ( ) 

{ 

top ( )  ; 

text ("FRAGMENT  SPRAY  ZONE"); 

input 2 ( ) ; 

input21() ; 

loc_2 ( ) ; 

vector2 ( ) ; 

vector2 1 ( ) ; 

index2 ( ) ; 

output 2 ( ) ; 

output 2 1 ( ) ; 

} 


/**  SUBROUTINE:  DRAW3()  **/ 


draw3 ( ) 

{ 

top ( ) ; 
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text ("BEST  POINT  for  DETONATION"); 

loc_3() ; 

vector3 () ; 

index3 ( ) ; 

output 3 () ; 

output 3 1 ( ) ; 

coord_3 ( ) ; 

repeat ( ) ; 

} 

/*********************************•****•****+***********•*/ 

/*     SUB-SUBROUTINE   INPUT1()         for  drawl ()        */ 

/it********************************************************/ 

input 1 ( ) 

{ 

varl4 () ; 

setcolor(&white) ; 

bxl=476 . 0 ;byl=165 . 0 ;bx2=638 . 0 ;by2=295 . 0 ; 

gcx=577. 0;gcy=294.0; 

boxf ( ) ; 

yarl3 () ; 

settextclr(&ligreen, Sblue) ; 

tcx=540. 0;tcy=282.0; 

movtcurabs ( &tcx, &tcy) ; 

text ("DATA") ; 

settextclr(&liyellow/&blue) ; 

tcx=510. 0;tcy=269. 0; 

movtcurabs (&tcx, &tcy) ; 

text ("Target") ; 

settextclr (Swhite, &blue) ; 

tcx=490. 0;tcy=258.0; 

movtcurabs (&tcx, &tcy) ; 

text("l.Tx  ="); 

tcy=247.0; 

movtcurabs(&tcx/ &tcy) ; 

text("2.Ty  =") ; 

tcy=23  6.0; 

movtcurabs ( &tcx , &tcy ) ; 

text("3.Vt  ="); 

tcy=225.0; 

movtcurabs ( &tcx , &tcy) ; 

text("4.Ht  ="); 

settextclr (Sliyellow, &blue) ; 

tcx=5 10 . 0 ; tcy=2 14.0; 

movtcurabs ( &tcx , &tcy ) ; 

text ("Missile") ; 

settextclr(&white,&blue) ; 

tcx=490. 0;tcy=203 .0; 

movtcurabs (&tcx, &tcy) ; 

text ("5. Mx  ="); 

tcy=192.0; 
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movtcurabs(&tcx,&tcy) ; 

text ("6. My  ="); 

tcy=181.0; 

movtcurabs (&tcx, &tcy) ; 

text("7.Vm  ="); 

tcy=170.0; 

movtcurabs ( &tcx, &tcy) ; 

text("8.theata=") ; 

} 

/*     SUB-SUBROUTINE    INPUT11()        for  drawl ()        */ 

/************************************** *******************/ 

inputll() 

{ 

settextclr ( Swhite , &blue) ; 

sprintf (s,"%-8. If", tx) ; 

tcx=555 . 0 ; tcy=258 . 0 ; 

varlO () ; 

sprintf  (s,"%-8.  If  'Sty)  ; 

tcy=247. 0; 

varlO () ; 

sprintf (s,"%-6. If", vt)  ; 

tcy=23  6. 0; 

varl0() ; 

sprintf  (s,n-5.1f»,ht)  ; 

tcy=225.0; 

varlO () ; 

sprintf (s,"%-8. If", mx)  ; 

tcy=203.0; 

varlO () ; 

sprintf (s,"%-8. If", my) ; 

tcy=192.0; 

varlO() ; 

sprintf (SjM%-6. If", vm)  ; 

tcy=181.0; 

varlO() ; 

sprintf (s,"%-5. If", theata) ; 

tcx=580 . 0 ; tcy=170 . 0 ; 

varlO() ; 

ht=ht*PI/180.0; 

theata=theata*PI/180 . 0 ; 

} 
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/*       SUB-SUBROUTINE    COORD_l()       for  drawl ()       */ 


coord_l ( ) 

{ 

gcx=60 . 0 ;gcy=60 . 0 ; 
movabs ( &gcx, &gcy) ; 
gcx=423 .5; 
lnabs(&gcx, &gcy) ; 
gcx=8  0.0; gcy=5  0.0; 
movabs (&gcx, &gcy) ; 
gcy=3  00. 0; 
lnabs(&gcx,&gcy) ; 

for  (count=0;count<6;count++) 

{ 

gcx=gcx+40*l. 439 ;gcy=62 . 5 ; 

movabs (&gcx, &gcy) ; 

gcy=60.01; 

lnabs(&gcx, &gcy) ; 

} 
gcy=60. 0; 
for  (count=0;count<6;count++) 

{ 

gcx=  8  0 . 0 1 ; gcy=gcy +  40.0; 

movabs ( &gcx , &gcy) ; 

gcx=84. 0; 

InabsC&gcx, &gcy) ; 

} 

settextclr(&liyellow) ; 
tcx=2  3  0.0; tcy=4 1.0; 
movtcurabs(&tcx/ &tcy) ; 
text("x  [ft]") ; 
tcx=3  0.0; tcy=2  00.0; 
movtcurabs (&tcx, &tcy) ; 
text("y  [ft]"); 
settextclr(&lired) ; 
tcx=350. 0;tcy=37.0; 
movtcurabs (&tcx, &tcy) ; 
text("[scale=200]")  ; 
} 
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/Ik********************************************************/ 

/*      SUB-SUBROUTINE    LOC_l()        for  drawl ()         */ 
/******************************************************** */ 

loc_l() 

{ 

var3 ( ) ; 
acftl() ; 
msl()  ; 
} 


ft********************************************************/ 

/*     SUB-SUBROUTINE   VECTOR1 ( )      for  draw()  */ 

ft********************************************************/ 

vectorl ( ) 

{ 

dist_factor=abs(sx)/ (cos(atan(abs(sy/sx) ) ) ) ; 
dist_l=dist_factor*trans_f actor/2 . 5 ; 
dist_2=dist_factor*trans_f actor; 

/***********  vector   for  Viti        **********/ 

setcolor(Stlicyan)  ; 

dx=dist_l*cos(theata) ;dy=dist_l*sin(theata) ; 

dx=dx*1.439; 

lnrel (&dx, &dy) ; 

/***********    vector  for  Vt   **********/ 
setcolor(&white) ; 
if  (  ht  ==  0.0  ) 

dx=-dist_l*vt/vm; 
else 

dx=dist_l*vt/vm ; 
dy=0 . 0 ; 
dx=dx*1.439; 
lnrel (&dx,&dy) ; 

/****•******    vector  for  Vmt  **********/ 

Inqgcur ( &gcx , &gcy , &white ) ; 

gcxl=gcx ; gey l=gcy ; 

setcolor(&liyellow) ; 

gcx=imx+80 . 0 ;gcy=imy+60 . 0 ; 

movabs(&gcx, &gcy) ; 

lnabs(&gcxl, &gcyl) ; 

/*********   extension  of  Vmt   **********/ 
if  (  ht  ==  0.0  ) 

dx=dist_2*cos(theata) -dist_2*vt/vm; 
else 


42 


dx=dist_2*cos (theata) +dist_2*vt/vm; 

dy=dist_2*sin(theata) ; 
dx=dx* 1.439; 
style=3 ; 

setlnstyle(&style) ; 
lnrel(&dx, &dy) ; 
style=l; 

setlnstyle(&style) ; 
} 


/*      SUB-SUBROUTINE    INDEX1 ( )      for  drawl ()        */ 
/••••••••A************************************************/ 

indexl ( ) 

{ 

dx=22. 0;dy=0.0; 
gcx=12  0. 0;gcy=3  00.0; 
tcx=150 . 0 ; tcy=296 . 0 ; 
in_vm() ; 

gcy=290.0; 
tcy=28  6. 0; 
in_vt ( ) ; 

gcy=2  80.0; 
tcy=27  6.0; 
in_vmt() ; 

gcy=270.0; 
tcy=2  66.0; 
in_miss() ; 
} 


/•••A****************************************************  */ 

/*      SUB-SUBROUTINE    0UTPUT1()       for  drawl ( )        */ 
/it********************************************************/ 

output 1 ( ) 

{ 

setcolor (Swhite) ; 

bxl=476 . 0 ;byl=70 . 0 ;bx2=638 . 0 ;by2=163 . 0 ; 

gcx=577 . 0 ;gcy=162 . 0 ; 

boxf ( ) ; 

settextclrf&ligreen, &blue) ; 
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tcx=535 . 0 ; tcy=152 . 0 ; 

movtcurabs (&tcx,&tcy) ; 

text ("OUTPUT") ; 

settextclr(&white,&blue)  ; 

tcx=490.0;tcy=141.0; 

movtcurabs (&tcx,&tcy) ; 

text ("9.  Sx   =") ; 

tcy=13  0.0; 

movtcurabs ( &tcx , &tcy ) ; 

text("10.Sy   ="); 

tcy=119.0; 

movtcurabs ( &tcx , &tcy ) ; 

text("ll.Vmtx  =") ; 

tcy=108.0; 

movtcurabs (&tcx, &tcy) ; 

text("12.Vmty  =") ; 

tcy=97.0; 

movtcurabs (&tcx,&tcy) ; 

text("13.tau  ="); 

tcy=86.0; 

movtcurabs ( &tcx, &tcy) ; 

text ( " 14 . minsx=" ) ; 

tcy=75.0; 

movtcurabs (&tcx, &tcy) ; 

text("15.minsy=") ; 

} 


/*     SUB-SUBROUTINE    0UTPUT11()     for  drawl ()  */ 

/*********************************************************/ 

output 11 ( ) 

{ 

sprintf (s,"%-9.1f",sx) ; 

tcx=57  2 . 0 ; tcy=14 1 . 0 ; 

varlO ( ) ; 

sprintf (s,"%-9. If", sy> ; 

tcy=130.0; 

varlO() ; 

sprintf (s ,"%-7. 2f"/vmtx) ; 

tcy=119.0; 

varlO() i 

sprintf (s,"%-7.2f"/vmty) ; 

tcy=108.0; 

varlO() ; 

sprintf (s,M %-6. 4 fM,tau) ; 

tcy=97.0; 

varlO() ; 

sprintf (s/"%-7.2f"/minsx) ; 

tcy=8  6.0; 
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varlO() ; 

sprintf (s,"%-7. 2 f",minsy) ; 

tcy=75.0; 

varlO ( ) ; 

sprintf (s , "%-7. 3f",m_dist) ; 

h=2;w=l; 

settext(&h,&w,&p,&m)  ; 

tcx=3  60 . 0 ; tcy=3  2  0.0; 

varlO ( ) ; 

} 


/•••A*****************************************************/ 

/*     SUB-SUBROUTINE   0UTPUT12()      for  drawl ()         */ 
/*********************************************************/ 

output 12 ( ) 

{ 

/******      miss  distance     **********/ 

setcolor(&liblue) ; 
gcx=itx+80 . 0 ;gcy=ity+60 . 0 ; 
movabs (&gcx, &gcy) ; 
dx=-minsx* 1.439  *trans_f actor ; 
dy=-minsy *trans_f actor ; 
lnrel (&dx, &dy) ; 

/**  determine  for  early  or  late  bird  ***/ 

settextclr(&green, &black) ; 

tcx=3  4  0.0; tcy=2  80.0; 

movtcurabs (&tcx, &tcy) ; 

deltcur() ; 

if ( ( (ty>my) &&( ( (ht==0. 0) &&(atan(vmty/vmtx) <atan(sy/sx) ) ) 

|  |"((ht>3.12)&&(atan(vmty/vmtx)>atan(sy/sx))) )) 

I  I 

( (ty<my) &&( ( (ht==0.0) && (atan(vmty/vmtx) >atan(sy/sx) ) ) 

I | ((ht>3.12)&&(atan(vmty/vmtx)<atan(sy/sx) ) ) ) ) ) 

text ("Early  Bird  !") ; 
else 

text ("Late  Bird  !") ; 
gwrite(filel) ; 
bottom () ; 
varl2 ( ) ; 
} 


/*    SUB-SUBROUTINE   INPUT2 ( )  for  draw2()       */ 

/•••••a************************************** 

input 2 ( ) 
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{ 

alphal=alphal*180 . O/PI ; 
alpha2=alpha2*180 . O/PI ; 

varl4 ( ) ; 

setcolor(&white) ; 

bxl=460 . 0 ;by 1=195 . 0 ;bx2=638 . 0 ;by2=305 . 0 ; 

gcx=577 . 0 ;gcy=299 . 0 ; 

boxf ( ) ; 

varl3 ( ) ; 

settextclr(&ligreen, Sblue) ; 
tcx=5  32.0; tcy=2  91.0; 
movtcurabs ( &tcx, &tcy) ; 
text ("INPUT") ; 

settextclr(&liyellow,&blue) ; 
tcx=4  95.0; tcy=2  80.0; 
movtcurabs ( &tcx , &tcy) ; 
text ("Encounter") ; 
settextclr(&white,&blue) ; 
tcx=4  70.0; tcy=2  70.0; 
movtcurabs (&tcx,&tcy) ; 
text("l.R_det  =") ; 
settextclr(&liyellow,&blue) ; 
tcx=495 . 0 ;tcy=260 . 0 ; 
movtcurabs (&tcx, &tcy) ; 
text ("Warhead") ; 
settextclr(&white,&blue) ; 
tcx=4  7  0.0; tcy=2  50.0; 
movtcurabs (&tcx, &tcy) ; 
text("2.alphal  =") ; 
tcy=240.0; 

movtcurabs (&tcx, &tcy) ; 
text("3.alpha2  =") ; 
tcy=23  0. 0; 

movtcurabs ( &tcx , &tcy) ; 
text ("4.N(f rag) =") ; 
tcy=22  0.0; 

movtcurabs ( &tcx , &tcy ) ; 
text ("5.  Vo    =") ; 
settextclr(&liyellow, &blue) ; 
tcx=495 . 0 ; tcy=210 . 0 ; 
movtcurabs (&tcx,&tcy) ; 
text ("Target") ; 
settextclr(&white, &blue) ; 
tcx=470.0;tcy=200.0; 
movtcurabs  (Stcx^tcy)  ; 
text ("6. length  =") ; 
} 
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/it********************************************************/ 

/*     SUB-SUBROUTINE     INPUT21()         for  draw2()      */ 
/•••••A***************************************************/ 

input21() 

{ 

sprintf (s,"%-5.1f",r_det) ; 

tcx=555 . 0 ; tcy=270 . 0 ; 

varlO() ; 

sprintf  (s,"%-6.  If  ",alphal)  ; 

tcy=250. 0; 

varlO () ; 

sprintf (s,"%-6. If", alpha2) ; 

tcy=240.0; 

varlO () ; 

sprintf  (s,"%-7. if" ,num)  ; 

tcy=2  3  0.0; 

varlO () ; 

sprintf (s,"%-7. If", vo) ; 

tcy=220.0; 

varlO () ; 

sprintf (s,"%-5. If", length) ; 

tcy=2  00.0; 

varlO() ; 

} 


/A*******************************************************/ 
/*    SUB-SUBROUTINE    LOC_2 ( )         for  draw2()        */ 

loc_2 ( ) 

{ 

var7 ( ) ; 
msl() ; 
} 


/A********************************************************/ 
/*    SUB-SUBROUTINE    VECTOR2 ( )       for  draw2()         */ 


vector2 () 

{ 

alphal=alphal*PI/180 . 0 ; 

alpha2=alpha2*PI/180.0; 
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gcx=imx+8  0.0; gcy=imy+6  0 . 0 ; 
movabs(&gcx,&gcy) ; 

/********   vector  for  Vo  due  to  alphal  *********/ 

setcolor(&ligreen) ; 

dx=dist_l* (vo/vm) *cos (theata+alphal) ; 

dy=dist_l* (vo/vm) *sin (theata+alphal) ; 

dx=dx* 1.439; 

lnrel (&dx, &dy) ; 

/***********       vector  for  Vm      ***********/ 
setcolor(&licyan) ; 
dx=dist_l*cos(theata) ; 
dy=dist_l*sin(theata) ; 
var9 ( ) ; 

/************      vector  for  Vi      ***********/ 
setcolor(&lipurple) ; 

dx=dist_l*l. 439* (cos (theata)+ (vo/vm) *cos (theata+alphal) ) ; 
dy=dist_l* (sin (theata)+( vo/vm) *sin( theata+alphal) ) ; 
lnrel (&dx, &dy) ; 

/************      vector  for  Vt      ***********/ 

setcolor(&white) ; 

dx=-dist_l*vt/vm ; 

dy=0 . 0 ; 

dx=dx* 1.439; 

lnrel (&dx, &dy) ; 

/*************     vector  for  Vft    ************/ 

setcolor(&lired) ; 

inqgcur(&gcx, &gcy,&lired) ; 

gcxl=gcx ; gey l=gcy ; 

gcx=imx+80 . 0 ;gcy=imy+60 . 0 ; 

lnabs(&gcx, &gcy) ; 

} 


/*     SUB-SUBROUTINE     VECT0R21()      for  draw2()        */ 

vector21() 

{ 

/**********   vector  for  Vo  due  to  alpha2   **********/ 

setcolor(&ligreen) ; 

dx=dist_l* (vo/vm) *cos (theata+alpha2) ; 

dy=dist_l* (vo/vm) *sin(theata+alpha2) ; 
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dx=dx*1.439; 
lnrel (&dx,&dy) ; 

/*************     vector  for  Vm    ****************/ 
setcolor(&licyan) ; 
dx=dist_l*cos(theata) ; 
dy=dist_l*sin(theata) ; 
var9 ( ) ; 

/***************   vector  for  Vi   *****************/ 
setcolor(Slipurple) ; 

dx=dist_l*1.439*(cos(theata)+(vo/vm) *cos(theata+alpha2) ) ; 
dy=dist_l*(sin(theata)+(vo/vm) *sin(theata+alpha2) ) ; 
lnrel (&dx,&dy) ; 

/****************    vector  for  Vt   ***************/ 

setcolor(&white) ; 

dx=-dist_l*vt/vm ; 

dy=0 . 0 ; 

dx=dx*1.439; 

lnrel (&dx, &dy) ; 

/****************    vector  for  vft  ***************/ 

setcolor(&lired) ; 

inqgcur (Sgcx, &gcy, Slired) ; 

gcx2=gcx ; gcy2=gcy ; 

gcx=imx+80 . 0 ;gcy=imy+60 . 0 ; 

lnabs(&gcx, &gcy) ; 

/****     triangle  for  fragment  spray  zone     *****/ 

movabs ( &gcx2 , &gcy2 ) ; 

lnabs(&gcxl, &gcyl) ; 

settextclr(&lired, Sblack) ; 

tcx=15. 0;tcy=37.0; 

movtcurabs (&tcx, &tcy) ; 

text ("Triangle  with  red  lines  shows  fragment  spray  zone") ; 

tcy=25. 0; 

movtcurabs (&tcx, &tcy) ; 

text ("with  respect  to  target!"); 

} 


/•••••••A*************************************************/ 

/*     SUB-SUBROUTINE     INDEX2 ( )         for  draw2 ( )       */ 
/•••••••••••A*********************************************/ 

index2 ( ) 

{ 

dx=22.0;dy=0.0; 
gcx=2  0.0; gcy=3  00 . 0 ; 
tcx=47 . 0 ;tcy=296 . 0 ; 
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in_vo ( ) ; 

gcy=290.0; 
tcy=286.0; 
in_vm() ; 

gcy=2  80.0; tcy=2  76.0; 
in_vi ( ) ; 

gcy=2  70.0; tcy=2  66.0; 
in_yt ( ) ; 

gcy=260 . 0 ; tcy=256 . 0 ; 
in_vf t ( ) ; 

} 


/it********************************************************/ 

/*     SUB-SUBROUTINE     OUTPUT2()         for  draw2()      */ 

output 2 ( ) 

{ 

setcolor(&white) ; 

bxl=460 . 0 ;byl=47 . 0 ;bx2=638 . 0 ;by2=173 . 0 ; 

gcx=577 . 0 ;gcy=172 . 0 ; 

boxf  ( )  ; 

varl3 ( ) ; 

settextclr(&ligreen, &blue) ; 

tcx=53  0 . 0 ;tcy=160 . 0 ; 

movtcurabs ( &tcx, &tcy) ; 

text ("OUTPUT") ; 

settextclr(&white,&blue)  ; 

tcx=470 . 0 ; tcy=150 . 0 ; 

movtcurabs (&tcx, &tcy) ; 

text("l,Vftxl  =") ; 

tcy=140. 0; 

movtcurabs ( &tcx , &tcy ) ; 

text("2.Vftyl  ="); 

tcy=13  0.0; 

movtcurabs (&tcx,&tcy) ; 

text("3.Vftx2  ="); 

tcy=120.0; 

movtcurabs (&tcx, &tcy) ; 

text("4.Vfty2   ="); 

tcy=110.0; 

movtcurabs (&tcx, &tcy) ; 

text("5.Vftl   ="); 

tcy=100.0; 

movtcurabs (&tcx, &tcy) ; 
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text("6.Vft2    ='•); 

tcy=90.0; 

movtcurabs (&tcx, &tcy) 

text("7.Gammal  =") ; 

tcy=8  0.0; 

movtcurabs (&tcx, &tcy) 

text ( " 8 . Gamma2  = " ) ; 

tcy=70.0; 

movtcurabs (&tcx, &tcy) 

text("9.Phione  =") ; 

tcy=60.0; 

movtcurabs (&tcx, &tcy) 

text("10.Phitwo=") ; 

tcy=50.0; 

movtcurabs  ( &tcx,  &tcy) 

text("ll.Rho   =») ; 

) 


/•••••••••••A*********************************************/ 

/*      SUB-SUBROUTINE       OUTPUT21()      for  draw2 ( )     */ 
/***********************************^ 


output2 1 ( ) 

{ 

gammal=gammal*180 . 0/PI ; 
gamma2=gainma2  *180 .  0/PI ; 
phione=phione*180 . 0/PI ; 
phitwo=phitwo*180 . 0/PI ; 


sprintf 

tcx=556 

varl0() 

sprintf 

tcy=140 

varl0() 

sprintf 

tcy=13  0 

varl0() 

sprintf 

tcy=120 

varlO() 

sprintf 

tcy=110 

varlO() 

sprintf 

tcy=100 

varlO() 

sprintf 


s,"%-8.2f",vftxl)  ; 
0;tcy=150.0; 

s,w%-8*2fw,vftyl) ; 

0; 

s,"%-8.2f",vftx2) ; 
0; 

s,"%-8.2f",vfty2) ; 

s,"%-8.2f",vftl) ; 
0; 

s,"%-8.2f",vft2) ; 
0; 

s,H%-6.1fM,gammal)  ; 
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tcy=90.0; 

varlO()  ; 

sprintf  (s,  "%-6.  If '  ,gamma2)  ; 

tcy=80.0; 

varlO ( ) ; 

sprintf (s, "%-6. If '  ,phione) ; 

tcy=7  O.O; 

varlO() ; 

sprintf (s,"%-6. If", phitwo) ; 

tcy=60.0; 

varlO () ; 

sprintf (s,"%-7. 4  f',rho) ; 

tcy=50.0; 

varlO() ; 

gwrite(file2) ; 
bottom () ; 
varl2 () ; 

} 


/it********************************************************/ 

/*      SUB-SUBROUTINE       LOC_3()       for  draw3()       */ 
/************************************  * / 

loc_3 ( ) 

{ 

trans_f actor=0 . 02*1.2; 

var8 ( ) ; 

acf tl ( ) ; 

} 


/*      SUB-SUBROUTINE      VECTOR3()       for  draw3 ( )      */ 

vector3 ( ) 

{ 

/****   vector  for  Vft  from  tail  section  of  aircraft  ****/ 

dx=0. 0;dy=0.0; 

movrel (&dx, &dy) ; 

dx=-trans_f actor *vftxl*l . 439 ; 

dy=-trans_f actor*vf ty 1 ; 

lnrel(&dx,&dy) ; 

dx=-dx ; dy=-dy ; 

movrel (&dx, &dy) ; 

dx=-trans  factor*vftx2*1.439; 
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dy=-trans_f actor*vf ty2 ; 
lnrel (&dx, &dy) ; 
dx=-dx ; dy=-dy ; 
movrel (&dx, &dy) ; 

/****   vector  for  Vft  from  pitot  tube  of  aircraft    ****/ 
dx=84. 0;dy=0.0; 
if  (ht==0.0) 

dx=dx ; 
else 

dx=-dx; 
movrel (&dx, &dy) ; 
dx=-trans_f actor*vf txl*l . 439 ; 
dy=-trans_f actor*vf ty 1 ; 
lnrel (&dx, &dy) ; 
dx=-dx ; dy=-dy ; 
movrel (&dx,&dy) ; 
dx=- trans_f actor * vf tx2 * 1 . 4  3  9 ; 
dy=-trans_f actor*vf ty2 ; 
lnrel (&dx, &dy) ; 
mdist() ; 
vmt_vector() ; 
msl_pt() ; 
} 


/*     SUB-SUBROUTINE     INDEX3 ( )        for  draw3 ( )       */ 
/*********************************************************/ 

index3 () 

{ 

varl3 () ; 

dx=22 . 0;dy=0.0; 

gcx=5  2  0.0; gcy=7  0.0; 

tcx=547. 0;tcy=66.0; 

in_vf t ( ) ; 

gcy=6  0 . 0 ; tcy=5  6 . 0 ; 
in_miss() ; 

gcy=5  0.0; tcy=4  6.0; 

in_vmt ( ) ; 

} 


/*      SUB-SUBROUTINE      OUTPUT3()        for  draw3()     */ 
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output 3 ( ) 

{ 

setcolor(&white) ; 

bxl=520 . 0 ;by 1=110 . 0 ;bx2=639 . 0 ;by2=200 . 0 ? 

gcx=577 . 0 ;gcy=199 . 0 ; 

boxf  ( )  ; 

settextclr(&ligreen,&blue) ? 
tcx=540 . 0 ; tcy=185 . 0 ; 
movtcurabs(&tcx,&tcy) ; 
text("Tgt  length"); 

tcy=165.0; 

movtcurabs(&tcx,&tcy) ; 
text ("Miss  Dist") ; 

tcy=145. 0; 

movtcurabs(&tcx, &tcy) ; 
text("Gammal") ; 

tcy=125.0; 

movtcurabs(&tcx, &tcy) ; 
text  ("Gamnia2")  ; 
} 


/*     SUB-SUBROUTINE     0UTPUT31()      for  draw3()      */ 
/••ye******************************************************/ 

output31() 

{ 

settextclrC&white^blue)  ; 

sprintf (s,"%-5. If" , length) ; 

tcx=560 . 0 ; tcy=175 . 0 ; 

varlO  (■)  ; 

sprintf (s ,"%-?. 3f"/m_dist) ; 

tcy=155.0; 

varlO () ; 

sprintf  (s,"%-6.  If",  gairraial)  ; 

tcy=135.0; 

varl0() ; 

sprintf (s,"%-6. If", gamma2) ; 

tcy=115.0; 

varl0() ; 

} 
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/it********************************************************/ 

/*  SUB-SUBROUTINE     COORD_3()         for  draw3()      */ 

/**•****************•*************************************/ 

coord_3 ( ) 

{ 

gcx=2  0.0; gcy=4  0.0; 
movabs(&gcx,&gcy) ; 
gcx=440.0; 
lnabs(&gcx,&gcy) ; 
gcx=2  0.0; gcy=4  0.0; 
movabs ( &gcx , &gcy) ; 
gcy=272. 0; 
lnabs(&gcx,&gcy) ; 

for  (count=0;count<5;count++) 

{ 

gcx=gcx+84 . 0 ;gcy=42 . 5 ; 
movabs (&gcx, &gcy) ; 
gcy=40.01; 
lnabs(&gcx, &gcy) ; 

} 
gcy=40.0; 
for  (count=0;count<4 ;count++) 

{ 

gcx=2 0.01; gcy=gcy+84/l. 439; 
.  movabs ( &gcx, &gcy) ; 
gcx=24.0; 
lnabs(&gcx, &gcy) ; 
} 

settextclr(&ligreen) ; 
tcx=4  50.0; tcy=3  8.0; 
movtcurabs ( &tcx, &tcy) ; 
text("x  [ft]") ; 
tcx=3 . 0 ; tcy=2  80.0; 
movtcurabs ( &tcx , &tcy ) ; 
text("y  [ft]") ; 
tcx=200.0;tcy=27.0; 
movtcurabs ( &tcx , &tcy ) ; 
text ( " [ scale=tgt  length] " ) ; 
deltcur() ; 
} 


/•••••A***************************************************/ 

/*     SUB-SUBROUTINE     REPEAT ( )         for  draw3()     */ 
/••••••A**************************************************/ 

repeat ( ) 


55 


{ 

gwrite(file3) ; 
iter_l : 
bottom () ; 
bottoml ( ) ; 
iter_3 : 
c=getch() ; 
if  (c=='f ) 

{ 

iter_2 : 

setcolor (&black) ; 
clr(); 

gread(file2) ; 
bottom () ; 
bottoml ( ) ; 
iter_4 : 
c=getch() ; 
if  (c=='f) 

{ 

setcolor(&black) ; 
clr(); 

gread(filel) ; 
bottom () ; 
c=getch() ; 
goto  iter_2 ; 

} 
else 

{ 

setcolor (Sblack) ; 
clr()  ; 

gread(file3) ; 
goto  iter_l; 
} 
} 
else 

{ 

setcolor (&black) ; 

clr(); 

close () ; 

} 


/*****         MISCELLANEOUS  SUBROUTINE  FUNCTION       *****/ 
/it********************************************************/ 


/**********         aircraft  function         **********/ 
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acftl() 

{ 

varl4 () ; 

setcolor(&white) ; 
gcx=itx+80 . 0 ;gcy=ity+60 . 0 ; 
movabs ( &gcx, &gcy) ; 
setcolor(&liyellow) ; 

dx=3.0;dy=0.0; 

var2 ( ) ; 

dx=-6 . 0 ; dy=0 . 0 ; 

var2 ( ) ; 

dx=3 .0;dy=0.0; 

var2 ( ) ; 

dx=0 . 0 ; dy=-3 . 0 ; 

var2 ( ) ; 

dx=0. 0;dy=6.0; 

var2 ( ) ; 

dx=0 . 0 ; dy=-3 . 0 ; 

var2 ( ) ; 

setcolor (Swhite) ; 

if  (ht==0.0) 

{ 

gcx=itx+3  6. 0;gcy=ity+56. 0; 

} 
else 

{ 
gcx=itx+l24.0;gcy=ity+56. 0; 

} 

movabs ( &gcx, &gcy) ; 

dx=0. 0;dy=3.0; 

var2 ( ) ; 

dx=2 . 0;dy=0.0; 

var2 ( ) ; 

dx=2.0;dy=16.0; 

var2 ( ) ; 

dx=6.0;dy=0.0; 

var2  ()  ; 

dx=4 . 0 ; dy=-12 . 0 ; 

var2 ( ) ; 

dx=12  .  0 ;  dy=-l .  0 ; 

var2 ( ) ; 

dx=14.0;dy=3.0; 

var2 ( ) ; 

dx=10.0;dy=4.0; 

var2 ( ) ; 

dx=4.0;dy=0.0; 

var2 ( ) ; 

dx=10 . 0 ; dy=-6 . 0 ; 

var2 ( ) ; 
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dy=-1.5; 
dy=-2 . 5 ; 
dy=0 . 0 ; 
dy=-0.5; 


dx=8 . 0 

var2 ( ) 

dx=8 . 0 

var2 ( ) 

dx=4 . 0 

var2 ( ) 

dx=0 . 0 

var2 ( ) 

dx=-4 . 0 ; dy=0 . 0 ; 

var2 ( ) ; 

dx=-8 . 0 ; dy=-2 . 5 ; 

var2 ( ) ; 

dx=-8 . 0 ; dy=-l . 0 ; 

var2 ( ) ; 

dx=-2  4 . 0 ; dy=0 . 0 ; 

var2 ( ) ; 

dx=-4 . 0 ; dy=-2 . 0 ; 

var2 ( ) ; 

dx=-6.0;dy=0.0; 

var2 ( ) ; 

dx=-10 . 0 ; dy=2 . 0 ; 

yar2 ( ) ; 

dx=-18 . 0 ; dy=l . 0 ; 

var2 ( ) ; 

dx=-2 . 0 ; dy=0 . 0 ; 

var2 ( ) ; 

setcolor(Slired) ; 

dx=-9 . 0 ;dy=0 . 25 ; 

var2 ( ) ; 

dx=9 . 0 ; 

var2 ( ) ; 

dx=-7 . 5 ; 

var2 ( ) ; 

dx=7 . 5 ; 

var2 ( ) ; 

dx=-8 . 4 ; 

var2 ( ) ; 

dx=8 . 4 ; 

var2 ( ) ; 

dx=-9.1; 

var2 ( ) ; 

dx=9 . 1 ;. 

var2 ( ) ; 

dx=-9 . 0 ; 

var2 ( ) ; 

dx=9 . 0 ; 

var2 ( ) ; 

} 

/************+  itiissile   function  ****************/ 

msl() 
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{ 

varl4 ( ) ; 

setcolor(&white) ; 
gcx=imx+80 . 0 ;gcy=imy+60 . 0 ; 
movabs (&gcx, &gcy) ; 
msl_point() ; 

gcx=imx+80 . 0 ;gcy=imy+60 . 0 ; 
movabs (&gcx, &gcy) ; 
msl_in_l_point() ; 
} 


msl_point() 

{ 

dx=10.0;dy=0.0; 
var6 () 


0 ; dy=2 . 0 ; 
0;dy=0.0; 
dy=-2 . 0 ; 


0;dy=-2.0; 
0;dy=0. 0; 
dy=2 . 0 ; 


dx=-10 

varl ( ) 

dx=-30 

varl() 

dx=0 . 0 

varl() 

dx=40.0;dy=0. 0; 

var6 () 

dx=-10 

varl() 

dx=-30 

varl() 

dx=0.0 

varl() 


setcolor(&lired) ; 

dx=-2  0 . 0 ; dy=0 . 0 ; 

varl() 

dx=2  0.0; dy=0 . 5 ; 

var6 () 

dx=-18 

varl() 

dx=18 . 0 ;dy=-l . 0 ; 


0;dy=0.0; 


var6() 

dx=-18 

varl() 

dx=18 . 0 ; dy=l . 5 ; 

var6 () 

dx=-16 

varl() 

dx=16 . 0 ;dy=-2 . 0 ; 

var6() 

dx=-16 

varl() 

dx=16 . 0 ; dy=2 . 5 ; 


0;dy=0.0; 


0;dy=0.0; 


0;dy=0.0; 


missile  point  function    *************/ 
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var6() ; 

dx=-14 . 0 ;dy=0 . 0 ; 

varl ( ) ; 

dx=14 . 0 ; dy=-3 . 0 ; 

var6 ( ) ; 

dx=-14 . 0 ; dy=0 . 0 ; 

varl  ( )  ; 

} 

/************  missile  initial  point  function  ***********/ 
msl_in_l_point ( ) 

{ 

setcolor(&cyan) ; 
dx=-3.0;dy=0.0; 
lnrel ( &dx, &dy) ; 
dx=6.0;dy=0. 0; 
lnrel (&dx,&dy) ; 
dx=-3 . 0 ; dy=0 . 0 ; 
lnrel (&dx, &dy) ; 
dx=0. 0;dy=-3 .0; 
lnrel ( &dx, &dy) ; 
dx=0.0;dy=6.0; 
lnrel (&dx,&dy) ; 
dx=0. 0;dy=-3 . 0; 
lnrel ( &dx, &dy) ; 
dx=50 . 0 ; dy=0 . 0 ; 
lnrel (&dx,&dy) ; 
dx=-50 . 0 ;dy=0 . 0 ; 
lnrel (&dx, &dy)  ;* 

inqgcur(&gcx,&gcy,&cyan)  ; 
gcxl=gcx ; gey l=gcy ; 

dx=50 . 0*cos (theata) *1. 439 ; 
dy=50 . 0*sin (theata) ; 
lnrel (&dx,&dy) ; 
movabs(&gcxl, Sgcyl) ; 

radius=2  0.0; 

angl=0 . 0 ; ang2=theata ; 

arc ( Sradius , Sangl , &ang2 ) ; 

movabsC&gcxlz&gcyl) ; 

} 


/***********  miss  distance  vector  function  *************/ 
mdist() 

{ 

setcolor(Stliblue)  ; 
gcx=itx+80 . 0 ; gcy=ity+60 . 0 ; 
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movabs (&gcx, &gcy) ; 

dx=-84 . 0*minsx/ length; 

dy=- (84 . 0*minsy/ length) /l. 439 ; 

lnrel (&dx,&dy) ; 

} 


/*********    missile  velocity  vector  function  **********/ 
vmt_vector ( ) 

{ 

setcolor(&liyellow) ; 
dx=vmtx*  2  * trans_f actor * 1.439; 
dy=vmty*2  *trans_f actor ; 
lnrel (&dx, &dy) ; 


dx=-2. 0*dx; 
dy=-2. 0*dy; 
lnrel (&dx,&dy) ; 

} 


/*********    variable  missile  point  function    *********/ 
msl_pt ( ) 

{ 

varl4 () ; 

setcolor(&white) ; 
inqgcur(&gcx,&gcy, &white) ; 
gcxl=gcx ; gey l=gcy ; 
msl_point () ; 
movabsC&gcxl, &gcyl) ; 
msl_in_l_point() ; 
} 


/*************         index  of  Vft         **************/ 
in_vf t ( ) 

( 

setcolor(&lired) ; 

movabs(&gcx, &gcy) ; 
lnrel (&dx, &dy) ; 
settextclr (&lired, &black) ; 
movtcurabsC&tcx, &tcy) ; 
text ("Vft") ; 
} 


/*************  index  of  vmt         *************/ 

in_vmt ( ) 

{ 

setcolor(&liyellow) ; 
movabsC&gcx, &gcy) ; 
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lnrel(&dx,&dy) ; 
settextclr (Sliyellow, Sblack) t 
movtcurabs (&tcx, &tcy) ; 
text("Vmt") ; 

} 


/**************    index  of  miss  distance    ***************/ 
in_miss ( ) 

{ 

setcolor(&liblue) ; 
movabs (&gcx, &gcy) ; 
lnrel(&dx, &dy) ; 
settextclr(&liblue, Sblack) ; 
movtcurabs (&tcx, &tcy) ; 
text ("Miss  Dist") ; 
} 


/**************         index  of  vi         **************/ 
in_vi ( ) 

{ 

setcolor(&lipurple) ; 

movabs ( &gcx, &gcy) ; 

lnrel (&dx, &dy) ; 

settextclr (&lipurple, &black) ; 

movtcurabs (&tcx, &tcy) ; 

text("Vi") ; 

} 


/************     index  of  missile  velocity   ************/ 
in_vm ( ) 

{ 

setcolor(&licyan) ; 

movabs (&gcx, &gcy) ; 

lnrel (&dx,&dy) ; 

settextclr (&licyan,&black) ; 

movtcurabs ( &tcx, &tcy) ; 

text ( "Vm" ) ; 

} 


/*************     index  of  target  velocity   ************/ 
in_vt ( ) 

{ 

setcolor(&white) ; 
movabs (&gcx, &gcy) ; 
lnrel (&dx,&dy) ; 
settextclr ( Swhite , &black) ; 
movtcurabs (&tcx, &tcy) ; 
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text("Vt") ; 
} 


/*************  index  of  fragment  velocity   ***********/ 

in_vo ( ) 

{ 

setcolor(&ligreen) ; 
movabs (&gcx, &gcy) ; 
lnrel(&dx/ &dy) ; 
settextclr(&ligreen, Sblack) ; 
movtcurabs ( &tcx, &tcy) ; 
text("Vo") ; 
} 

/***************  top  marking        **************/ 

top() 

{ 

h=2 ;w=2 ; 

settext ( &h , &w , &p , &m) ; 

settextclr (Slicyan, &black) ; 

tcx=2  0. 0; 

tcy=320.0; 

movtcurabs ( &tcx , &tcy) ; 

} 


/***************      bottom  marking      ***********.*****/ 
bottom ( ) 

{ 

width  =  13; 

setlnwidth(&width) ; 

setcolor(&liblue) ; 

gcx=10 . 0 , gcy=10 .  0 ; 

movabs(&gcx, &gcy) ; 

gcx=63  5.0; 

lnabs(&gcx, &gcy) ; 

varl3() ; 

settextclr (Swhite, &liblue) ; 

tcx=450. 0;tcy=7. 0; 

movtcurabs (&tcx, &tcy) ; 

deltcur() ; 

text ("Hit  any  key  to  proceed") ; 

> 


bottoml ( ) 

{ 

varl3 () ; 

settextclr(&white,&liblue) ; 
tcx=150 . 0 ; tcy=7 . 0 ; 
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movtcurabs ( &tcx , &tcy ) ; 

text ("Hit") ; 

tcx=198.0; 

movtcurabs ( &tcx, &tcy) ; 

text ("to  see  previous  display.")? 

settextclr(&lired,&liblue) ; 

tcx=182.0; 

movtcurabs (&tcx, &tcy) ; 

deltcur() ; 

text("F") ; 

} 


/****************     routine  box  function  **************/ 
boxf ( ) 

{ 

box ( &bxl , &byl , &bx2 , &by2 ) ; 
movabs(&gcx, &gcy) ; 
flood (&blue) ; 
} 


/****************     closing  screen    ******************/ 
close () 

{ 

h=7 ; w=5 ;p=0 ;m=0 ; 

settext (&h,&w, &p, &m) ; 

settextclr (Slicyan, Sblack) ; 

tcx=240. 0;tcy=14  0. 0; 

movtcurabs ( &tcx, &tcy) ; 

deltcur() ; 

text ("END") ; 

varl2 ( ) ; 

} 


/***  routine  varing  functions  ***/ 

/it********************************************************/ 

varl() 

{ 

if (((3.14<theata)&&(theata<4.71) ) | | ( (1.57<theata) && (theata< 
3.14))) 

{ 

dx=dx*cos(theata) -dy*sin (theata) ; 

dy=dx* sin (theata) +dy*cos (theata) ; 

dy=-dy ; 

Inrel (&dx, &dy) ; 
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} 

else 

if (((0.0<theata)&&(theata<1.57) ) | | ( (4 .71<theata) && (theata<6 

.28))) 

{ 

dx=dx*cos (theata) -dy*sin (theata) ; 

dy=dx*sin (theata) +dy*cos (theata) ; 

lnrel (&dx, &dy) ; 

} 

} 

var2 ( ) 

{ 

if  (ht==0.0) 

dx=dx ; 
else 

dx=-dx; 

lnrel (&dx, &dy) ; 
} 

var3 () 

{ 

trans_f actor  =0.2  ; 
if ( (txonx) && (ty<my) ) 

{ 

itx=0 . 0 ; ity=0 . 0 ; var5 ( ) ; 

} 
else  if ( (tx<mx) && (ty>my) ) 

{ 

itx=0.0;ity=2  00. 0;var5 () ; 

} 
else  if ( (tx>mx) && (ty>my) ) 

{ 

imx=0 . 0 ; imy=0 . 0 ; var4 ( ) ; 

} 
else  if ( (tx>mx) &&(ty<my) ) 

{ 

imx=0. 0;imy=200.0;var4 () ; 

} 
} 

var4 ( ) 

{ 

itx=imx+ (tx-mx) *trans_f actor* 1 . 439 ; 

ity=imy+ (ty-my) *trans_f actor ; 
} 

var5() 

{ 

imx=itx+ (mx-tx) *trans_f actor*l . 439 ; 
imy=ity+(my-ty) *trans_f actor; 
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} 

var6 ( ) 

{ 

if (((3.14<theata)&&(theata<4.71) ) | | ( (1.57<theata) && (theata< 

3.14))) 

{ 

dx=dx*cos (theata) -dy*sin (theata) ; 

dy=dx* sin (theata) +dy*cos (theata) ; 

dy=-dy; 

movrel(&dx, &dy) ; 

} 

else 

if ( ( (0.0<theata) &&(theata<1.57) ) | | ( (4.71<theata) &&(theata<6 

.28))) 

{ 

dx=dx*cos (theata) -dy* sin (theata) ; 

dy=dx* sin (theata) +dy*cos (theata) ; 

movrel (&dx, &dy) ; 

} 

} 


var7 ( ) 

{ 

if  (((tx<mx)&&(ty<itiy) )  &&  (alphaKO.  0) ) 

{ 

imx=2  80.0; imy=4  0.0; 

} 
else  if  ( ( (tx<mx) &&(ty<my) ) && (alphal>0. 0) ) 

{ 

imx=120 . 0 ; imy=180 . 0 ; 

} 
else  if  (( (tx<mx)  &&(ty>my)  )&&  (alphaKO.  0) ) 

{ 

imx=0 . 0 ; imy=10 . 0 ; 

} 
else  if  ( ( (tx<mx) &&(ty>my) ) && (alphal>0. 0) ) 

{ 

imx=280 . 0 ; imy=110 . 0 ; 

} 
else  if  (( (tx>mx)  &&(ty>my)  )&&  (alphaKO.  0) ) 

{ 

imx=10 . 0 ; imy=160 . 0 ; 

} 
else  if  ( ( (tx>mx) &&(ty>my) ) && (alphal>0.0) ) 

{ 

imx=140 . 0 ; imy=10 . 0 ; 

} 
else  if  (( (tx>mx)  &&(ty<my)  )&&  (alphaKO. 0) ) 

{ 
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imx=2  00.0; imy=15  0 . 0 ; 

} 

else  if  ( ( (tx>mx) &&(ty<my) ) && (alphal>0. 0) ) 

{ 

imx=4  0.0; imy=4  0.0; 

} 
} 

var8 ( ) 

{ 

if  (((tx<mx)&&(ty<my) ) &&(alphal<0.0) ) 

{ 

itx=12  0.0; ity=13  0 . 0 ; 

} 
else  if  ( ( (tx<mx) &&(ty<my) ) && (alphal>0. 0) ) 

{ 

itx=13  0 . 0 ; ity=2  0.0; 

} 
else  if  ( ( (tx<mx)  &&(ty>my) )  &&  (alphaKO.  0) ) 

{ 

itx=19  0 . 0 ; ity-17  0.0; 

} 
else  if  ( ( (tx<mx) &&(ty>my) ) &&(alphal>0.0) ) 

{ 

itx=12  0 . 0 ; ity=8  0.0; 

} 
else  if  (( (tx>mx)  &&(ty>my)  )&&  (alphaKO.  0) ) 

{ 

itx=2  6  0 .  0  ;  ity=7  0 .  0 ; 

} 
else  if  ( ( (tx>mx) &&(ty>my) ) && (alphal>0. 0) ) 

{ 

itx=22  0 . 0 ; ity=150 . 0 ; 

} 
else  if  (( (tx>mx)  &&(ty<my)  )&&  (alphaKO.  0) ) 

{ 

itx=110 . 0 ; ity=3  0.0; 

} 
else  if  ( ( (tx>mx) &&(ty<my) ) && (alphal>0. 0) ) 

{ 

itx=110 . 0 ; ity=160 . 0 ; 

} 

} 

var9 ( ) 

{ 

dx=dx* 1.439 ; 

lnrel (&dx, &dy) ; 

gcx=imx+80 . 0 ;gcy=imy+60 . 0 ; 

movabs(&gcx, &gcy) ; 

} 


67 


varlO() 

{ 

movtcurabs (&tcx, &tcy) ; 

text(s) ; 

deltcur() ; 

} 

varl2() 

{ 

c=getch() ; 
setcolor(Stblack)  ; 
clr(); 

} 


varl3() 

{ 

h=l ;w=l ;p=0 ;m=0 ; 

settext ( &h , &w , &p , &m) ; 

} 

varl4() 

{ 

width=l ; 

setlnwidth(Swidth) ; 
} 


varl5() 

{ 

if  (  atan(sy/sx)  <  atari  (vmty/vmtx)  ) 

{ 

alphal  =  -alphal  ; 

alpha2  =  -alpha2  ; 

} 
else 

{ 

alphal  =  alphal  ; 
alpha2  =  alpha2  ; 
} 


varl6() 

{ 

if  ( ( (vftxKO.O)  &&(vftyl>0.0)  )  |  |  (  (vftxKO.O)  &&(vftyl<0.0)  )  ) 

gamma l=gamma 1+3 . 14 ; 
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else 

gamma l=gamma 1+6 .  28 ; 
if  (gammal>6.28) 

gamma l=gammal-6 . 28 ; 
else 

gamma l=gammal ; 

if  (((vftx2<0.0)&&(vfty2>0.0)) | | ((vftx2<0.0)&&(vfty2<0.0))) 

gamma2=gamma2+3 . 14 ; 
else 

gamma2=gamma2+6 . 28 ; 
if  (gamma2>6.28) 

gamma2=gamma2-6 . 28 ; 
else 

gamma2=gamma2 ; 
} 
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APPENDIX  C 


USER'S  MANUAL 


The  following  is  the  user's  manual  of  ENDGAME 
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USER'S  MANUAL 


AE  3251/AE  3705 
AIR  DEFENSE  LETHALITY 


A  STUDY 

of 

the  ENDGAME 

in 
AIR  DEFENSE 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 
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I .    INTRODUCTION 

This  homework  problem  is  given  to  the  student  to 
present  an  opportunity  for  him/her  to  see  how  the 
vulnerability  of  an  aircraft  can  be  assessed  in  a  missile- 
based  air  defense  situation.  The  computer  simulation 
program  named  ENDGAME  will  be  used  to  assess  the  missile 
miss  distance,  the  fragment  spray  zone  and  the  best  point 
for  detonation. 


II.   PROBLEM  DEFINITION 

You  are  going  to  conduct  an  endgame  assessment  of  a 
generic  attack  aircraft  on  an  attack  mission  to  destroy  a 
ship  or  ground  vehicle.  Each  student  will  be  given  an 
ordnance  package,  target  and  missile  velocity  and  flight 
path.  This  problem  is  purely  for  instructional  purposes 
and  is  not  based  on  any  actual  or  planned  combat  situation. 
The  aircraft  vulnerability,  missile  characteristics, 
ordnance  package,  flight  path  parameter  limits  have  been 
chosen  only  to  provide  guidelines  for  the  class  problem. 


1.  The  initial  positions  of  the  target  and  the  missile 
are  restricted  to  within  100,000  ft  and  must  be  above 
ground.  That  is,  the  initial  x  and  y  positions  of  the 
target  and  the  missile  must  be  realistic.  The  error 
message  which  tells  you  "  UNREAL!  ASSUME  THAT  IT  BE  ABOVE 
GROUND  AND  REAL."  comes  if  you  type  incorrect  or 
unrealistic  input  data.   Just  input  reasonable  data. 


2.  The  relative  difference  in  distance  between  the 
initial  x  and  y  positions  of  the  target  and  the  missile 
must  each  lie  within  1,000  ft.  If  this  requirement  is  not 
met,  no  warning  message  will  appear,  but  the  results  on  the 
screen  may  be  meaningless. 


3 .  The  target  or  aircraft  speed  must  be  greater  than  0 
and  less  than  2000  ft/sec  (  or  Mach  2  )  .  Otherwise,  an 
error  message  tells  you  "  UNREAL!  " .  Just  type  again  the 
appropriate  input. 


4.  Assume  straight  and  level  flight  for  the  air 
target.  The  target  heading,  therefore,  must  be  either  0° 
(left  to  right)  or  180°  (right  to  left) . 
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5.  The  missile  speed  must  be  greater  than  the  target 
speed.  If  not,  the  error  message  "  UNREAL!"  appears.  So, 
use  reasonable  input  data. 

6.  The  input  for  the  elevation  angle  of  the  missile, 
6,  must  also  be  reasonable.  If  the  x  and  y  positions,  for 
example,  of  the  target  is  greater  than  those  of  the 
missile,  and  if  e  is  given  as  145°  or  210°,  then  the  error 
message  will  warn  you  to  use  the  appropriate  elevation 
angle  in  order  to  allow  the  target  to  be  approached. 


III.  EXECUTING  ENDGAME 

You  can  execute  ENDGAME  on  any  IBM  PC/AT  personal 
computer  or  compatible  with  a  monitor  graphic  resolution 
640  x  350.  If  you  have  a  different  system  on  a  personal 
computer,  a  few  modifications  to  the  program  are  required 
in  order  to  run  ENDGAME.  Further  details  describing  these 
modifications  are  given  in  this  user's  manual.  The 
following  instructions  tell  you  how  to  run  ENDGAME.  Three 
different  displays  will  come  out  from  ENDGAME  program. 


1.  Turn  on  the  system  power  of  your  personal  computer. 

2 .  Insert  the  diskette  which  contains  the  execution 
file,  the  source  file  and  the  auxiliary  files. 

3.  Run  the  executable  file,  endgame.exe,  by  simply 
typing  "ENDGAME".  Then,  you  will  see  the  opening 
screen. 

4.  Just  follow  the  instructions  shown  on  the  screen  to 
do  the  simulation. 

5.  Do  not  hit  any  key  to  proceed  until  the  message 
line  comes  on  the  bottom  of  the  each  screen.  This 
means  that  the  current  display  is  being  written  to 
the  specified  file  while  the  In-Use  light  on  the 
diskette  driver  comes  on. 

6.  You  will  be  given  the  opportunity  to  see  the 
previous  screen  by  hitting  "F"  key  when  the  best 
point  for  detonation  screen  is  displayed. 
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IV.   OUTPUT 

Make  sure  that  your  printer  is  ready  if  you  want  to 
obtain  the  hard  copy  of  the  display.  Just  hit  the  "prt  sc" 
key  while  holding  the  shift  key,  or  hit  the  "print  screen" 
key  by  itself  on  other  special  computer  keyboards. 


V.    WAYS  TO  MODIFY  THE  PROGRAM 

To  modify  the  program  in  order  to  use  the  system 
that  you  have  or  you  are  going  to  use,  you  must  make  the 
following  changes. 

1.  Among  the  external  variables  declared,  the  number 
of  mode  and  page  which  are  declared  by  integer  should  be 
varied  by  the  systems  that  you  have.  That  number  depends 
upon  how  much  memory  is  installed  on  the  board  and  what 
kind  of  graphic  device  driver  you  have.  See  section  7  in 
Function  Description  of  HALO  manual  (Version  2.26). 

2.  The  graphics  device  driver  must  also  be  changed. 
See  the  HALO  manual. 
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